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In recent measurements of tuBneling transport through individual ferromagnetic Co nanograins, 
Deshmukh, Gueron, Ralph et alSia (DGR) observed a tunneling spectrum with discrete resonances, 
whose spacing was much smaller jfehan what one would expect from naive independent-electron esti- 
mates. In a previous publication, El we had suggested that this was a consequence of nonequilibrium 
excitations, and had proposed a "minimal model" for ferromagnetism in nanograins with a discrete 
excitation spectrum as a framework for analyzing the experimental data. In the present paper, we 
provide a detailed analysis of the properties of this model: We delineate which many-body electron 
states must be considered when constructing the tunneling spectrum, discuss various nonequilib- 
rium scenarios and compare their results with the experimental data of Refs. We show that 
a combination of nonequilibrium spin- and single-particle excitations can account for most of the 
observed features, in particular the abundance of resonances, the resonance spacing and the absence 
of Zeeman splitting. 

PACS numbers: 73.23.Hk, 75.50.Cc, 73.40.Gk 



I. INTRODUCTION 

An important milestone in the study of itinerant fer- 
romagnetism was reached during the last two years, 
when Deshmukh, Gueron, Ralph et alxM (DGR), us- 
ing single-electron tunneling spectroscopyjj succeeded for 
the first time to resolve discrete resonances in the tun- 
neling spectrum through individual ferromagnetic single- 
domain Cobalt nanograins, with diameters between 1 and 
4 nm. Their work goes beyond previflus-|Sl;udies of fer- 
romagnetic single electron transistorsj3l3'LlEl which eluci- 
dated the interplay of ferromagnetism and charging ef- 
fects: the fact that DGR's Co grains were sufficiently 
small that discrete resonances could be resolved, means 
that they were probing the true quantum states partici- 
pating in electron tunneling, which allows the nature of 
electron correlations in itinerant ferromagnets to be stud- 
ied in unprecedented detail. Besides the intrinsic scien- 
tific interest of studying ferromagnetism on the nm-scale, 
the insights so gained might also be of technological inter- 
est, since the size of memory elements in magnetic storage 
technologies is decreasing extremely rapidlyjj and papti- 
cles as small as 4 nm are coming under investigation.E3 

An examination of the magnetic-field dependence of 
the individual resonances observed by DGR indicated 
that two-current models, in which spin- up and spin-down 
electron bands are considered effectively independently, 
are inadequate for describing the true electronic states 
inside a nm-scale ferromagnet - so that a fundamen- 
tally different theoretical approach is required. To this 
end, a simple new phenomenological model was recently 
introduced by thejjaresent authors together with Desh- 
mukh and Ra,lph,B and independently by Canali and 
MacDonald.EJ We regard this as a "minimal model" for 
ferromagnetic nanograins, in that it seems to be the sim- 
plest model possible for a discrete-state system which 
takes into account the electronic correlations induced by 



magnetic interactions and which treats the particle's tc^ 
tal spin as a quantum mechanical variable. We arguedtl 
that the main features of the measured tunneling spectra 
can be understood within this model by assuming that 
nonequilibrium spin accumulation occurs, which can pro- 
duce a much denser spectrum of tunneling excitations 
than expected within an independent-electron model. 
The conclusion that nonequilibrium effects play an im- 
portant role has since|-|been confirmed by more recent 
measurements by DGRa on a gated device, in which new 
resonances appeared as the gate voltage was tuned to 
drive the system further away from equilibrium. 

In the present paper, we provide a detailed analysis of 
equilibrium and nonequilibrium tunneling through ferro- 
magnetic grains, within the framework of our minimal 
model. The present analysis goes well beyond that of 
Ref. H in that we consider not only the spin ground 
states of Fig. 2 in Rcf. |^, but also spin wave excitations 
and single-particle excitations. The latter turn out to be 
necessary to understand why a small resonance spacing 
is observed even when the threshold bias voltage for the 
onset of tunneling is rather small, so that nonequilibrium 
effects cannot be very strong. 

The structure of the paper is as follows. In Sec. II 
wej-£ummarize the expjerimental results of Deshmukh et 
al.a and Gueron et alB The Hamiltonian of our minimal 
model is presented in Sec. together with a convenient 
set of basis states for analyzing the low-lying excitations 
and their energies. In Sec. IV we discuss two different 
equilibrium excitations in ferromagnetic grains, namely 
single particle excitations and spin excitations. A de- 
tailed discussion of nonequilibrium excitations and their 
consequences for tunneling spectra is given in Sec. 0: 
First we calculate the current through a grain for a 
nonequilibrium scenario involving only transitions be- 
tween the ground states \s,s) of a ladder of spin mul- 
tiplets of different total spin s; the resonance spacing for 



2 



the peaks in the conductance is found to be in quantita- 
tive agreement with DGR's measurements if we assume 
a total ground state spin of about sq ~ 1000 and that 
the resonances are predominantly due to the tunneling 
of minority electrons. We then generalize this nonequi- 
librium scenario by including also all more highly-excited 
states \s,m) of these spin multiplets, finding that if the 
total spin is large (s ^ 1), the Zeeman splitting of the 
observed resonances should be strongly suppressed, in 
agreement with DGR's experiments. Finally, we show 
that a combination of single-particle- and spin excita- 
tions in the presence of nonequilibrium can lead to the 
large number of resonances seen in experiments. Some 



concluding remarks can be found in Sec. VI. 

A brief account of our results on the most simple 
nonequilibrium scenario in ferromagnetic grains has al- 
ready been given in Ref. ^. The present pape r includes a 
detailed derivation of these results, (Sec. VB ), since this 
paves the way for the mor e co mplicate d non equilibrium 
scenarios presented in Sec. V_C and Sec. VP , which have 
not been reported before. 



II. SUMMARY OF EXPERIMENTAL RESULTS 

In DGR's experiments,0i a nanoscale Cobalt grain was 
used as a central island in a single electron transistor: it 
was connected via tunnel barriers to external leads and 
for one of the grains in Ref. ^, the central grain was also 
capacitively coupled to a gate. The electronic spectrum 
of the particle was determined by measuring the tunnel 
conductance through the grain as function of transport 
voltage {V), [gate voltage (V^)] and magnetic field fioH 
at a fixed temperature of < 90mK. The diameters of the 
Co-grains were estimated as 1-4 nm. Assuming a roughly 
hemispherical shape, the number of atoms in each grain 
then is in the range iVa « 20-1500, implying a total spin 
of So « 0.83iVa « 17-1250. 

Since the charging energy (> 30mcV) was very much 
larger than typical values of the transport voltage {eV < 
9meV) and the temperature, fluctuations in electron 
number on the grain are strongly suppressed, so that co- 
herent superposition between different electron numbers 
A'' need not be considered. The energy balance condition 
that determines through which eigenstates of the grain 
electrons can tunnel for given values of transport (and 
gate) voltage thus involve differences between eigenenep-. 
gies of a grain with fixed particle number A'' or ± 1 ,tHl 
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FIG. 1: The spin wave multiplet |s,m)Q for s = 1, A = 6. 
The position of each level represent its kinetic energy Sj. 



The excitation spectra measured by DGR had several 
properties that differ strikingly from those-|0f previously- 
studied non-magnetic Al and Au graina-H13 including: 
(PI): Many more low- energy excitations were observed 
than expected: For all Co grains studied, the ob- 
served level spacing is dobs S 0-2 meV, which is much 
smaller than the independent-electron estimate of d^in ~ 
1.2eV/so|-aBd dmaj ~ 4.6eV/so for minority-/majority 
electrons, EJtj respectively. (In Ref. ^ many more reso- 
nances were observed than in Ref. |l|, due to lower-noise 
tunnel barriers.) 

(P2): In the small-field regime {fioH < 0.2 T), the tun- 
neling resonance energies show strong non-linear depen- 
dence on H, with hysteresis and jumps at the switching 
field \Hsw\- For fields beyond the switching field, the lev- 
els can exhibit non-monotonic variations as a function of 
H, with each level behaving differently. For a more de- 
tailed discussion of the magnetic field dependence of the 
tunneling energies see Refs. |,| 

(P3): In the large-field regime {\H\ 3> |iJsw|), the reso- 
nances depend roughly linearly on H, with 7?-slopes that 
almost all have the same sign for a given grain, i. e. slopes 
of opposite signa,|iue to Zeeman-shifting of spin-up and 
spin-down levelsoEj are not observed. 
(P4): Measurements on a gated devicecl showed that 
the observed resonances correspond predominantly to the 
tunneling of minority electrons. p. 
(P5): Aleasuremcnts on a gated devicen also showed that 
some of the resonances must be due to nonequilibrium 
excitations, since some resonances disappear when the 
Coulomb blockade threshold for the onset of tunneling 
(i.e. the amount of nonequilibrium) is reduced by tuning 
the gate voltage to lie close to a degeneracy point. 



AE%^Ef^'-Ef, (1) 

each corresponding to the energy cost of some rate- 
limiting electron tunneling process \i)^ — ^ |/)^^^ onto 
or off the grain. Here \i)^ denotes a discrete eigenstate, 
with eigenenergy , of a grain with A^ electrons, etc. 
As the magnetic field H is swept, the resonances undergo 
energy shifts and crossings. 



III. MODEL 



In this section we introduce a simple model for ultra- 
small ferromagnetic grains. The challenge is to describe 
the individual quantized electronic excitations of a fer- 
romagnetic nanoparticle, taking into account the elec- 
tronic correlations induced by magnetic interactions and 
anisotropy forces. 
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Hamiltonian 



We propose to model a nanoscale magnet with disojete 
excitations by the following "minimal" HamiltonianJj'EEl 
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He + Hq + HcKch + Hzcc + Huni , (2a) 
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eVoSN + Ec5N 
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with h — QcsfJ-Bl^oH ■ Here He is the standard Coulomb 
charging energy for a nanoparticle with SN excess elec- 
trons. Vd = [Qo + ErVrCr] /C, for r = {L,R,g), 
is the electrostatic potential energy of the island, with 
C = J2r^r- Cli Cu and Cg are the tunnel junctions 
connecting the grain to the left or right lead or the gate 
electrode, respectively, and Vl, Vr and Vg are the volt- 
ages of the left/right leads and the gate electrode. Qo is 
an initial random off-set charge and Eq = e^/2C is the 
unit of charging energy. 

Hq describes the kinetic energy of a single band of 
single-electron states \j, a) , labeled by a discrete index j 
and a spin index a = (f , i), with the spin quantization 
axis chosen in the z-direction. The exchange, Zeeman, 
and anisotropy terms, Tioxch, Hzcc and Tiuni, are func- 
tions of the total spin vector S . 
spin vector of the electrons in level j, 



X] , Sj , where Sj is the 



(3) 



((T are Pauli matrices). Ticxch is a rotationally invari- 
ant term which models the effects of an exchange field 
and forces the system to adopt a non-zero total ground 
state spin, say sq. On account of this term, spins 
aligned parallel or antiparallel to {S) may be thought 
of as forming "majority" and "minority" bands, which 
effectively rotate rigidly together with the magnetiza- 
tion direction. We shall take the mean level spacings 
near the respective Fermi energies, dmin ~ 1.19eV/so 
and c?,„aj ~ 4.61 eV/so, and the exchange-splitting of t 
Fermi energies, Ap = ep.maj — £F,min (~ 2 eV for Co), 
as characteristic parameters of the model. 

Hicc describes the spin Zeeman energy in an exter- 
nal magnetic field H = Hz. Finally, the uniaxial 
anisotropy TYuni is the simplest non-trivial form of an 
anisotropy modeling the combined 
shape and surface anisotropics, etcEI 
tor in the easy-axis direction and fcAr(> 0) is a volume- 
independent constant, can be estimated from the 
measured switching field using fcjv ~ IJi-oIJi-bHsw, which 
yields Un ~ O.OlmeV [see Ref. 14 in Ref. |^. For com- 
pleteness we note that our analysis in a previous paper 



'ects of crystalline, 
fi is the unit vec- 



[Ref. H, see also Ref. || showed that the anisotropy con- 
stant fcjv undergoes fluctuations, i.e. varies between dif- 
ferent electronic states within a ferromagnetic nanopar- 
ticle. It was shown that these fluctuations have signifi- 
cant consequences for the magnetic field dependence of 
the resonance energies. However, in the present paper 
we shall focus only on understanding the characteristic 
spacings between resonances observed in DGR's experi- 
ments, and not on their magnetic field dependence; we 
will therefore neglect fluctuations in throughout this 
paper. 



B. Basis states 

Let us for the moment set the anisotropy strength 
fcjv = 0, and use the eigenstates of H{kN — 0) to con- 
struct a convenient set of "bare" basis states. Since the 
Hamiltonian then commutes with the total spin, these 
states can be grouped into spin multiplets that are la- 
beled by their S ■ S and S^ eigenvalues, say s{s + 1) and 
m. For example, the bare ground state of H{kN = 0) for 
given N, s and h (> 0), say \s, s)^ , can be written explic- 
itly as a member of the following multiplet of normalized 



states, \s,m)^ [illustrated in Fig. |l|: 
\s,s)^ = l[c]^ J|c]^ I vac) 



(4) 



s,to)o = 



Y (2,s)!(s-m)! ^'^-^ l'''^o ' 



Here — N/2 ± s, and 5_ = ^ ■ c^jCj^ is the spin- 
lowering operator. Within this model, the energy dif- 
ference between the spin multiplet |s,to) and all other 
states that are constructed from the same single-electron 
levels is at least otprder ep maj — min , a very large value 
(~ 2 eV for Co).^ 

The inclusion of a non-zero anisotropy term, Huni, will 
cause the true low-energy eigenstates, \s,m)^ , to be lin- 
ear superpositions of the bare states in the multiplet 



\s,m)!^ . We choose labels such that |s,m)^ 



as kN/h — > O"*". We shall call the states |s,to)^ the spin 
wave multiplet, since each can be viewed as a homoge- 
neous spin wave. 



|s,m)o 



C. Eigenenergies 

In the absence of anisotropics {Huni = 0), it is possi- 
ble to write down explicit expressions for the low-lying 
excitation energies of our model. Let us denote the ex- 
citation energy of a state |s,m)^ relative to the ground 
state |so,wo)o'" (for a gated device) byEi 



SE{6s,6m,SN) + eVD6N 



(6) 
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with Ss = s — So, Sm = m — toq and SN = N — Nq. It contribution to Eq. has the following form, first writ- 
is straightforward to show that the voltage-independent ten down by Canali and MacDonald:Eil 



SE{Ss,SmJN) = (SN)'^ 
+ds 



^maj ^iniii 



Er 



5N 



£F 



Af - U{so + 1/2) + ^ „ ^ 



4 4 
f (5s SN 



f(5s)2 

+ [dmaj + dmin " U] 



'^maj ^min 



h6m 



(7) 



Here we introduced the average of majority- and 
minority-band Fermi energies e-p = (eF.maj +eF,min)/2. 
The stability of the ground state spin sq, [*-jS. the re- 
quirement that 5E{Ss,5m,SN) > 0, impliesO the rela- 
tion Af = U{so + l/2) + do, where do(~ 1/so) is a small, 
grain-dependent energy satisfying 



A. Single particle excitations 



(dn 



U/2) + h<dQ < dn^in - U/2 + h 



(8) 



Hence, the magnitude of U may be estimated as J7 ~ 
Af/so — 2eV/so. Note that the ground state spin sq 
can be changed by a sufficiently large change in ap- 
plied magnetic field. However, the range over which 
the applied magnetic field has to be swept between 
two successive changes in ground state spin is of order 
Sh ~ (dmaj + dinin — U)/2, which is large in nanomagnets 
{fioSH > 25 T for Co particles with diameters < 4 nm). 
Therefore, for a given value of N, we shall, as long as we 
neglect nonequilibrium effects, consider only the "ground 
state" spin value sq. 



Apart 
Sec. fHB 



from the multiplets |s. 



discussed in 



higher energy multiplets can be built by cre- 
ating additional single-particle excitations, e.g. by start- 
ing from the bare multiplet constructed by applying the 
spin-lowering operator to the state c^ji-^c^ ^ |s, s)^, with 
'^i < j 5: '^T f ^ "T- However, their eigenenergies 
lie higher than those of the spin wave multiplet |s,to)^ 
by an amount of order Sji — Ej] this is at least of order 
the single-electron level spacing, dmaj — 4.61 eV/so (re- 
spectively c?inin = 1.19eV/so), i-c. rather large compared 
to c?obs [cf. (PI)]; thus the mechanism causing the ob- 
served abundance of low-energy excitations, whatever it 
is, cannot involve only purely single-particle excitations, 
but most also involve spin excitations. 



B. Spin wave excitations 



IV. EQUILIBRIUM TRANSITIONS 

To construct the tunneling spectrum associated with, 
say, adding an electron to the grain, we must in principle 
calculate the excitation energies ^E^- for all the allowed 
transitions, i.e. those for which the tunneling matrix el- 
ement 

Mj:^^+\f\cl\^f (9) 

is nonzero. For now, we shall neglect nonequilibrium ef- 
fects and thus consider only those tunneling processes for 
which the initial state corresponds to the grain's ground 
state, i.e. = |so,so)^. In particular, we shall fo- 
cus on two different types of equilibrium transitions: 
(A) transitions involving single-particle excitations whose 
resonance spacing, estimated from the electron density 
of states, is found to be much larger than observed in 
DGR's experiments and (B) transitions between differ- 
ent spin wave states, of which only two transitions are 
found to have significant weight, leaving unexplained the 
large number of resonances seen in experiments. 



We shall now study transitions between two multiplets 
of spin wave states, with initial states |z)^ — \si,mi)^ 
and final states I/) — \sf,mf) . Consider first the 
large-field regime h > /igw, [where /isw = ffoffM-BMo^^sw]- 
Since here Tizcc dominates over Tiuni, we may set = 
and construct the matrix elements Alj^ using the bare 

spin wave multiplets \si,mi)^ and \sf,mf)^^^. The 
condition Mj"^ ^ then implies the spin selection rules 
|s/ — Si| = 1/2 and \mf — mi\ = 1/2. 

Among all possible final states |/) satisfying these 
selection rules. Table | lists those three which can be 
reached by adding an electron to the lowest available lev- 
els of I So, so)q , namely ji = + 1 and j2 = ni + 1: A 
spint-electron can be added only to level ji (Table ||, row 
1), whereas a J,-electron can be added to either level ji 
(row 2) or j2 (row 3). 

The excitation energies AEj. (Table |, column 4) of 
the (ji t) and (ji |) transitions are degenerate at ft, = 
and Zeeman-split as a function of h, but, in accord with 
(P3), this splitting will not be observable: the weight of 
the (ji I) transition is smaller than that of the (ji t) 
transition (Table |, column 3) by a Clebsch-Gordan co- 
efficient of order l/(2so) which is negligibly small for 
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TABLE I: Matrix elements Mj^ and excitation energies AiJ^^ of those final states 1/)'^^^ = \sf,mf)^^^ that can be reached 
from the initial state = \so, so)^ by adding a spin-a electron to level ji = rij + 1 or j2 = n| + 1. 





Sf,mf 






AE+ 


iji T) 


so + 1/2, so + 1/2 


1 


Ec + EF.maj + dmaj 


^{U/2){so+3/4)^h/2 


Ui I) 


so + 1/2, so -1/2 


1/V2S0 + 1 


Ec + EF.maj + dmaj 


-([//2)(so + 3/4)+V2 


O2 i) 


so - 1/2, so - 1/2 


1 


Ec + ep.min + dmin 


+ (f//2)(so + l/4) + V2 



large-so grains.S The (ji t) and (j2 i) transitions both 
have large, comparable weights, and would produce reso- 
nances with large-h slopes of opposite signs. Depending 
on whether the difference in their excitation energies (Ta- 
ble |, column 4) is close to or far from (it is at most of 
order dimj — U/2, i.e. < 3meV for a 4-nm -diameter Co 
particletfl), either both or only one of the (ji t) and (j2 i) 
transitions would be observable in the regime of low- 
est excitation energies (say < O.SmeV). However, in an 
equilibrium tunneling scenario, this leaves unexplained 
the large observed density of tunneling resonances (PI), 
since, apart from these two transitions, there are no oth- 
ers with significant weight and excitation energies less 

than dmin [dme^j)- 

Next we argue that this problem persists also for the 
case with anisotropy, where ^ 0. A priori one might 
have expected to see more low-lying excitations then, 
since the selection rule \mf —mi\ = 1/2 no longer applies 
if the matrix elements Mj"^ are formed using the exact 

spin wave multiplets \so,mi)''^ and |so ± 1/2, m/)^'''^, 
which are not eigenstates. We have therefore numeri- 
cally diagonalized 7i zcc + '^uni Ji3 as function of h/k^, for 
a few selected values of sq, to determine the eigenstates 
\so,mi)^ and |so ± 1/2, m/)^'*'^ and calculate the ma- 
trix elements Afjf [Fig. |^, for both spin-increasing and 
-decreasing transitions (s/ = sq ± 1/2, j ~ "^1/1 + -'^)- 
We find that in both cases, the transition probability 

\M%\^ from |so, sq)"^ to |so ± 1/2, ± 1/2 - 
is very much larger for n = than for any other n ^ 
state. This is the same trend as that found in Table |. 
Thus, even though Tiuni causes violations of one of the 
spin selection rules, the extra transitions have too little 
weight to explain the large [density of low-energy excita- 
tions that is observed (Pl)£^ 

Apart from the fact that only two of the above dis- 
cussed transitions have significant weight there are two 
additional important considerations which lead us to con- 
clude that the abundance of resonances seen in experi- 
ments cannot be explained by equilibrium spin wave tran- 
sitions alone: (i) Firstly, the resonances associated with 
final states |s/,m/) that differ only in rtif would have a 
spacing of order fc^r [~ 0.01 meV], i.e. much smaller than 
the observed resonance spacing, (ii) Secondly, for high 
magnetic fields these resonances would exhibit a system- 
atic increase in the magnitude of their slopes (which are 
oc |si — TO/I) that was not observed in experiment. We 
therefore assert that the large density of resonances can- 
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FIG. 2: Matrix elements for the transitions |si,Si)^ 

.^Tfij)^^^ (si = 20), calculated for Tiuni + tLzcc and 
plotted as functions of /i/fcjv sweeping positive to negative, 
for the following final states: (a): |si — |,Si — |) (b): 



Si + 



1\JV+1 



(solid lines) and I Si + i , Si 



(dot- 



The solid lines are thick for M' 



ted lines). Compared to these, all other final states have neg- 
ligible matrix elements. 

hi 

\^)^\f))■ 



and thin for M^f (both contribute to the same transition 



not be explained by equilibrium transitions alone; we will 
explore nonequilibrium effects below. 



V. NONEQUILIBRIUM TRANSITIONS 

Since the large density of resonances (PI) cannot be 
explained by equilibriu m tra nsitions [neither single par- 
ticle excita tions (Sec. [V A ) nor spin wave excitations 
(Sec. IVB)], we shall in this section explore nonequilib- 
rium effects. Our conclusion will be that a combination of 
nonequilibrium spin and single-particle excitations pro- 
duces a much denser spectrum of tunneling states than 
expected within an independent-electron model. For a 
spin of So ~ 1000, this nonequilibrium scenario gives a 
resonance spacing of ~ 0.2meV for the spacing of res- 
onances due to the tunneling of minority electrons, in 
accord with the observed resonance spacing. 

Nonequilibrium spin accumulation had of course al- 
ready been studied previously in the context of single^ 
electron transistors with ferromagnetic componentsj3'E£l 
and spin accumulation for nanograins with disctete en- 
ergy levels was first analysed by Barnas et alE3 How- 
ever, these analyses all employed a single-particle descrip- 
tion in which all states that were considered were simple 
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FIG. 3: Illustration of a nonequilibrium scenario involving 
only spin ground states, ||5s;(5Ar) = |so + 5s^so + Ss)'^"'^^'^ , 
for the case that the first electron that tunnels is a minority 
electron (a = —1) that enters the grain {p — +1). Here 5s 
and 5N characterize the spin with respect to the overall spin 
gro und s tate: 5s = s — so and 5N = N — No (see Sees. Ill C 
and 
AK 



VB) 



The vertical arrows indicate the energy differences 
„ J , AiS"'^, and A_Eo°* . Charging transitions are numbered 



as n = 0, 1, 2 • 



Slater determinants of single-particle states. Within our 
present model, we have to go beyond this simple pic- 
ture by considering the true many-body eigenstates of 
the Hamiltonian, which are in general linear combina- 
tions of Slater determinants. 

After explaining the general idea of noneq uilibrium 
processes in ferromagnetic grains in Sec. VA, we shall 
describe different nonequilibrium scenarios. For each, 
we calculate the corresponding theoretical tunnel spectra 
and compare resonance spacings and the number of reso- 
nances with DGR's measurements. We show in Sec. VB 



that nonequilibrium spin excitations lead to resonance 



spacings as observed in measurements and in Sec. VC 
that Zeeman splittings are suppressed for large spin sq. 
A co mbin ation of spin- and single particle excitations 
(Sec. VD) significantly enhances the number of tunnel- 
ing resonances achievable for a given setting of the gate 
voltage, making it possible to explain the large number of 
resonances observed by DGR even for the case of a small 
Coulomb blockade threshold {i.e. weak nonequilibrium). 



sition from the iV-electron ground state |7g)^ to the 
{N ± l)-electron ground state |7g)^^^, it may be ener- 
getically possible for the subsequent "discharging" tun- 
neling transition to return the particle to an excited A^- 
electron state |7o}^ instead of |7g}^, provided the ap- 
plied voltage is sufficiently large, eV > — . Like- 
wise, further charging and discharging transitions may 
allow any of a large ensemble of states to be occupied 
at higher and higher levels of an energy ladder, ter- 
minating only when an energy-increasing transition re- 
quires more energy than the applied voltage provides. 
As the voltage is increased, the total current (or con- 
ductance) may increase stepwise (or show peaks) when 
thresholds are crossed to allow higher-energy transitions 
up the nonequilibrium ladder, thereby changing the oc- 
cupation probabilities of the ensemble of nonequilibrium 
states and opening new tunneling channels. 

Let {7} be the set of all states involved in a nonequi- 
librium ladder of excitations, i.e. the set of all discrete 
states with a nonzero occupation probability ^(7) for a 
given bias voltage, gate voltage, temperature and mag- 
netic field. To find the occupation probability ^(7) for 
all the states I7) of the ladder one has to solve a normal- 
ization condition . P{^') = 1 and a stationary master 
equation of the formll3 



= 



E 



{s^yP(7') - Ey^P(7) 



(10) 



for each 7. The first (second) term in Eq. (|T^ describes 
the rate at which the probability of a given configuration 
increases (decreases) due to electrons tunneling onto or 
off the grain, the remaining terms are associated with 
electronic relaxation and spin-fiip relaxation on the grain, 
respectively. 

T,j-yi is the total tunneling-induced transition rate from 
initial state I7') to final state I7). Considering sequential 
tunneling only, it has the form 



E E ' 

■■=L,Rp=± 



(11) 



where S''"'', (S'' ,) involves the coherent transfer of an 
electron onto (from) the grain from (onto) lead r and is 
given by 



77 



r'^El(7l4l7')P/(i?7-^7'-eK), 



A. General master equation 



77' 



\{l\c,aH)\^f{E^-Ey+eVr) . 



In general, A^-electron states other than the ground 
state can be populated during the process of current 
flow, and|— this may affect the experimental tunneling 
spectrum.ll3E3 A simple scenario is illustrated in Fig. 0. 
Even if a first tunneling event causes a "charging" tran- 



Here f{E) — l/(e^/''^-^+l) is the Fermi function and eVr 
is the electrostatic potential energy difference between 
lead r and the grain, 



eVr = eVr — eVo 



(12) 
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The third and fourth terms in Eq. (10) describe elec- 
tronic (spin-conserving) relaxation processes inside the 
grain. For simplicity, we shall only consider electronic re- 
laxation between energetically "neighboring" single par- 
ticle levels, i.e. we take 



= r' 



''El 



(7h 



It') ? 



(13) 



(Above, the notation \\5s;5N) is used.) The indices p 
and a are used to distinguish whether the first electron 
that tunnels enters or leaves the grain, p = (+1,-1), 
and whether it is a majority or minority electron, a — 
(maj,min) =1+1, -!)• 

Using Eq. (0), the threshold energy coats for transitions 



L"^i and L"^2 can be calculated to be:l 



(Generalizations of this assumption are straightforward, 
though cumbersome.) The last two terms of Eq. ( [lo| ) de- 
scribe the rate at which the probability of a given distri- 
bution increases (decreases) due to spin-flip relaxation in 
the ferromagnetic grain. For simplicity we shall assume 
all spin-flip relaxation rates to be much smaller than 
all other rates, T''*' < W,r<^\ and hence take S"^ = with 



throughout this paperS) (Again, it is straightforward 
to consider generalizations of this case.) Moreover, all 
rates F are assumed to be independent of the specific 
single particle level ia involved. 

|Xhe current trough the grain can then be calculated 
asB 



77' 



77 



77 



(14) 



B. Spin accumulation 

In a ferromagnetic particle, in addition to the nonequi- 
librium occupation of single-electroe, states discussed pre- 
viously for non-magnetic particles ,E3 nonequilibrium spin 
excitations are possible, too, if thej-spin-flip rate T^^ is 
smaller than the tunneling rate r*"".t3 In this case a lad- 
der of transitions will occur between states with different 
total spin s, causing each to have a finite occupation 
probaiji 1 i.ty . and thus leading to spin accumulation on the 
grain 

The simplest nontrivial case, namely a ladder of spin 
multiplet ground states Is, s) (see Fig. ||) was already 
briefly discussed in Ref. H. Below we shall discuss this 
case in more detail. We shall calculate the resonance 
spacing of steps in the current (or of peaks in the dif- 
ferential conductance) and the number of resonances for 
transitions between spin ground states. We shall charac- 
terize spin ground states by their spin and charge relative 
to the overall ground state, i.e. we shall write 



|5s;5iV) = |so + (Js, So + fc) 



No+SN 



(15) 

with 5s and 5N as defined in Sec. [II C after Eq. ||. We 
shall find that the resonance spacing agrees very well with 
DGR's measurements. 

Consider a sequence of nonequilibrium transitions 
forming a "ladder" (L"p) with "rung index" n — 
0,1,2,..., where each rung corresponds to a "charging" 
transition {L"\) followed by a "discharging" transition" 

(Lzy- 

\\apn;0) \\ap{n + 1/2)- p) ^ \ \ap{n + 1);0) . 



E{ap{n + 1/2)- p) - E{apn;Q) (16a) 

+ n [d„ - C//2] , 
E{ap(n + 1); 0) - E(ap{n + l/2);p) (16b) 
^EZ2 + [da - C//2] , 



AEZi - d„ (1 +p)/2 + .Be +p(eF + ado/2) -C//8, 
da(l - p)/2 -Ec- p{eF - ado/2) - 3 C//8 



AE"P 



(Above, the notation a means maj — min and min = 
maj.) Note that the total energy cost for the combined 
transitions Lq^ and 21 namely 



^ AEZ+AE-}, 
= d. 2 



— H da — h pado ^ 



(17a) 



is always > 0; this follows intuitively from the fact 
that AEZot is the excitation energy between the over- 
all ground state \\Ss;5N) = ||0;0) and the adjacent-spin 
ground state ||ap;0), and more formally from condition 
(|) on do. 

Assuming that the peaks in the conductance are due 
to successive charging transitions becoming accessible as 
the bias voltage is increased, the resonance spacin g fo r 
the ladder L°'P can readily be calculated using Eqs. (|l^): 



6EZ ^ AEZ+,,,-AEZ (18a) 
= d„ - U/2 . (18b) 

This result, which evidently depends only on whether the 
charging transition involves the tunneling of a majority 
or minority electron, a = (maj, min), can be intuitively 
understood as follows: The resonance spacing, SE"^, is 
a difference of energy differences, i.e. a type of (discrete) 
second derivative of the total energy. The contribution 
da refiects the discrete second derivative w.r.t. the quasi- 
particle number of the energy involved in creating, using 
only a-electrons, n particle-like or hole-like excitations 
relative to the overall ground state, and the term —U/2 
refiects the discrete second derivative w.r.t. the spin of 
the exchange energy. The partial cancellation between 
da and ~U/2 in Eq. (18b) reflects the opposite signs of 
the kinetic and exchange energies in the Hamiltonian (H), 
and is thus very generic. Very signiflcantly, since U/2 and 
dmin (but not d,„aj ) happen to be almost equal, this par- 
tial cancellation of two rather large energies produces a 
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much smaller energy, 

6E™ ~ (4.6 - 1.0)eV/so = 3.6 eV/so , (19a) 
SE^^'P - (1.2 - 1.0)eV/so = 0.2 eV/so , (19b) 

where 



100 



150 



200 



250 



and [II 



used the parameter estimates given in Sees. II 
^ This brings us to one of the main conclusions 
of this paper: the small resonance spacing of 0.2 m eV 
observed by DGR is consistent with prediction ( |l9h| ) for 
minority- electron charging events if the ground state spin 
is assumed to be about sq ~ 1000, which is within the 
estimated size range of DGR's grains. Satisfactorily, the 
conclusion that minority electrons dominate the charging 
transitions, which was. reached independently by Canali 
and MacDonald toojlj has recently been confirmed ex- 
perimentally for DRG's gated device, as has been the 
conclusion that nonequilibrium physics is involved [cf. 
points (P4) and (P5) of Sec. II]. 

We shall therefore henceforth consider only the case 
that charging transitions are due to minority electrons, 
i.e. we take a = min (but for notational brevity will 
sometimes still use the index a instead of "min" ) . The 
conductance will then show a limited number, say n^Jl,, 
of resonances due to charging transitions, with a rather 
small spacing of dmin ^ U/2, followed by an unlimited 
number of resonances due to discharging transitions, with 
a much larger spacing of dmaj — U/2. This can be seen 
as follows: let us consider for definiteness a circuit with 
Vl = —Vr — y/2, and suppose that eV > 0, so that elec- 
trons flow from left to right through the ferromagnetic 
grain. Then transition involves tunneling across the 
left junction ifp= 1, or the right junction if p = —1, [and 
likewise L"^2 involves tunneling across junction {R, L) for 
p = (—1, 1)]. The voltage thresholds, say eV^\ (or eV"^), 
needed to overcome the energy cost Ai?"*^ (or AE'"^) in 
order for the charging transition L"^^ (or the discharging 
transition L^^j) to occur, are determined by the condi- 
tions: 



peV, > AEZ 
peVp > AEZ 



(20a) 
(20b) 



Here eVp , the electrostatic potential energy difference be- 
tween lead p and the grain, is [from Eq. (^2|)] related to 
the actual applied voltage by 



eVp = p eVBp — eVg 



(21) 



where Bp = {Cp + Cg/2)/C is a capacitance ratio which 
converts applied voltage to energy,t3 and eVg = e{Qo + 
V^Ca)/C is an offset energy. It follows from Eqs. ( |20| ) 
and (|2l] ) that and ^'"^ given by 



BpeVZ = AEZ+peVg, 



peVg 



BpeVZ = Ai?3 



(22a) 
(22b) 



In particular, the threshold voltage V^^ for the first 
charging transition, Lq\^ determines the measured size 
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100 150 200 250 

eVBL(meV) 

FIG. 4: Current as function of cVBl for the following pa- 
rameters: a = min, p — +1, so = 10, dmin — 1.19eV/so, 
dmaj = 4.61eV/so, U = 2eV/so, Bl = 0.3, Br = 0.7, 
^thresh ^ 0.|^V/so, EZot = OeV/so, r^/r« = 0.8, = 0, 
r = 80 mK.c3 (a) Only spin excitations (b) Spin excitations 
and single-particle excitations: /Y^ = 10^. Arrows mark 
additional current steps due to combined spin and single- 
particle excitations. No significance should be attached to 
step heights here, since they depend on (unknown) tunnel, 
ing matrix elements, which for simplicity we took all equalEj 
Parameters are chosen according to Fig. 1 of Ref. 0. 



of the Coulomb-blockade region (in energy units), say 

^thresh . 



thresh 



eVo"fBp = AE^P+peV, 



g ' 



(23) 



This enables us to eliminat e p eV g from Eq. (|2^) and write 
it as [using Eqs. ( p^ ) and (17a)] 



Bp eVZ = n{d^ - U/2) + S*?-^'^ 



(24a) 



n{d^ - U/2) - i?^'"-*^ + AKf„, . (24b) 



Now, whenever £;t,hrcsh 



> 



1 A 



^tiircsh small, the inequality 



i.e. whenever 



(25) 



will hold for sufficiently small values of n (at least for 
n = 0). However, for large enough n it will cease to 
hold, since for the case a = min that we are considering, 
the "step size" dmaj — U/2 for eV^2 much larger than 
the "step size" d^i^ — U/2 for eV^\. Thus, the first few 
measured conductance resonances will be due to a se- 
quence of (rather closely-spaced) charging transitions, as 
opposed to (much more widely-spaced) discharging tran- 
sitions, because each time the bias voltage is incremented 
by eV"\ to make the next charging transition L'^^ ener- 
getically accessible, this bias voltage increment is already 
large enough (namely > el^"^) to also allow the discharg- 
ing transition L"^2 to occur. However, once the inequality 
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Eq. (^9) is violated, the subsequent discharging transi- 
tion Lj'^ will become possible only after the total bias 
voltage increment reaches eT^f , i.e. henceforth discharg- 
ing (instead of charging) transitions will determine the 
conductance resonances, which will henceforth be spaced 
much more widely. 

To calculate the total number of closely-spaced reso- 
nances due to charging transitions, njr^g, we must thus 
determine how large n can become before the condition 
Eq. ( p5| ) ceases to hold. Using Eqs. ( p^ and Bp + Bp — 1, 
this condition can be rearranged to yield an expression 
for n.jij'g, which is found to be given by the smallest integer 
larger or equal to 



';. |s„+1,s„-1> 
\ |s„+1,s„> 
\ |s„+1,s„+1 ) , 
n=1 



N.+1 

|s„+1/2!s„-3/2) 

|S„+1/2,S„-1/2)| 

|S„+1/2,S„+1/2)]| 



|S„,S„-2) 

|s„,s„-1 > 

|So,So> 



n=0 



K+1 



|s„-1/2,S-5/2> 
/|s„-1/2,S„-3/2>\ 
/|S„-1/2,S„-1/2>V 



|s„-1,^-3> / 

, |s„-1,^-2) / 

Is„-1,^-1> / 
n=1 



ap=1 



ap=-1 



^thresh 



Bp{d„ 



U/2) ~ Bp{d^ 



U/2) 



(26) 



The prediction that n'^^^ increases linearly with J^q"'^'''' 
is in qualitative agreement with Fig. 2 of Ref. How- 
ever, it is not quite consistent with more recent data 
on Co-grains,0 where, even when the Coulomb blockade 
threshold was very small (i?^'^''"'^'^ < 1 meV), the dif- 
ferential conductance showed many (> 10) peaks, i.e. 
many more than Eq. (|26| ) would predict. To illustrate 
this, we have solved Eq. ( p^ ) numerically for the dis- 
cussed transitions and calculated the current for the pa- 
rameters of the model and a threshold charging energy of 
^thresh ^ o.8eV/so as function of cVBl [see Fig. |(a)]. 
The calculated current shows three steps in the plotted 
region. Experimental data (see Fig. 1 in Ref. H) on the 
other hand show « 12 resonances for the parameters 
used. 

To summarize the conclusions of this section: a 
nonequilibrium scenario involving only spin ground 
states, and assuming charging transitions involving mi- 
nority electrons, results in a spacing of resonances, 
SEi-cs ~ 0.2 meV for sq = 1000, that agrees roughly with 
the spacing observed. However, the number of resonances 
predicted by Eq. ( p^ is sometimes much smaller than 
observed, namely when is small (on the order of 



C. Spin wave excitation 



FIG. 5: Illustration of a nonequilibrium scenario involving 
spin- wave states: Different multiplets \s, m) for No and A'o-l-1 
electrons on the grain are depicted. Due to selection rules only 
specific transitions (marked by arrows) are possible. 
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FIG. 6: Current as function of eVBh for the following pa- 
rameters: a — min, p = -fl, sq — 10, dmin = 1.19eV/so, 
d^aj = 4.61eV/so, U = 2eV/so, Bl = 0.4, Br = 0.6, 
kN = 0.01 meV, F^/F^ = 0.8, F^' = 0, F'=' = and a 
Coulomb blockade region of 7eV/so; (a-c): ft = meV; (d-f): 
h — 0.05 meV; (b,e) shows the fourth current step of (a,d) 
for T = 80mK, (c,f) shows the same step for a lower temper- 
ature of r = 0.8 mK. No significance should be attached to 
step heights here, since they depend on (unknown) tunneling 
matrix elements. 



In the preceding section, we considered only transi- 
tions between different spin ground states, |s, s); we shall 
now extend this scenario to include all higher lying states 
of the corresponding multiplets \s,m) (see Fig. g). We 
shall show that this results in a fine structure for the 
current steps which would be resolvable only for very 
low temperaturcaHa and hence would not be expected to 
be observable in DGR's present measurements. More- 
over, we find that the Zeeman splitting of resonances is 
strongly suppressed for a large spin sq ^ 1- 

Fig. ^ illustrates the transitions between different spin 
wave states of "neighboring" multiplets. For simplicity, 
we assume that both the magnetic field and the easy axis 



lie in the z-direction, so that the selection rule 1??^/ — 
rriil ~ 1/2 holds. 

We solved Eq. (|o|) numerically for a spinel of sq ~ 10 
and a Coulomb blockade region of 7eV/so [as in Fig. 2, 
sample 3 of Ref. 0| , and calculated the current for zero 
magnetic field [Fig. ^(a-c)], and non-zero applied field 
[Fig. I(d-f)]. 

Fig. ^(a) shows the current as function of energy 
(BlcV) for a temperature of T = 80 mK and zero 
magnetic field, h = eV. The current displays seven 
equally spaced steps. These steps belong to tran- 
sitions between successive sets of pairs of multiplets. 
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FIG. 7: (a) A Zeeman-split current step, where A/dc 
is the overaU height of steps due to all S^-decreasing transi- 
tions and A/ = A/docreasing + AJincroasing is the Overall current 
step, (b) AJdccrcasing/AJ, Calculated numerically for various 



so-values for the second step in the current and F /F 



0.1. 



e.g. the first one belongs to transitions between \sQ,m) 
and |so — 1/2, to'), the next one to transitions between 
|so — 1,to) and |so — 3/2, to'), etc. Fig. ||(b) shows the 
fourth current step of (a) on a finer energy scale. In 
Fig. ^(c) the same step is shown for a lower temper- 
ature of r = 0.8 mK, at which it now reveals sub- 
structure in the form of 14 finer steps. These small 
steps correspond to transitions between various states of 
the two multiplets, jsg — 3, m) and |so — 7/2, m), namely 
the first small step in (c) corresponds to a transition, 
|so-3,±(so-4)) ^ |so-7/2,±(so-7/2)), the next 
one to |so — 3, ±(so — 5)) \so — 7/2, ±(so — 9/2)), etc. 
The last step in (c) corresponds to jsp ~ 3, ±(so — 3)) 
So — 7/2, ±(so — 7/2)). It can be checked easily, us- 
ing Eq. (^), that their spacing is given by fcAr/so [— 
O.OlmeV/so]- Similarly, all other steps of (a) (except 
for the first one) have a substructure of smaller steps be- 
longing to all transitions between neighboring multiplets 
that are allowed by the selection rules (as indicated in 
Fig. ^). Note that the number of ( large -scale) steps in 
(a) is still given by Eq. (^ ) of Sec. VB , with a spacing 
given by dmin -U/2 0.2eV/so]- 

In Figs, md) to (f) a similar set of plots is shown 
as in Figs.^(a) to (c), but now in the presence of 
an applied magnetic field, h = 0.05 meV. Fig. ^(d) 
shows the current itself. Figs. 1(e) and (f) again show 
the fourth current step for two different temperatures. 
Fig. ^(e) shows that the step of (b) has Zeeman-split 
into two steps and (f) shows that these two steps cor- 
respond to two groups of transitions, namely all S*^- 
increasing transitions (steps on the left side) and S*^- 
decreasing transitions (right side) between the two mul- 
tiplets |so — 3,to) and |so — 7/2,to). The fact that (f) 
shows many more steps than (c) results from the fact 
that the applied magnetic field lifts the degeneracy of 



states |s, ±to). The first step in (f) belongs to a transition 
So — 3, So — 4) — > |so — 7/2, So — 7/2) the next one to 
So — 3, So ~ 5) — > |so — 7/2, so — 9/2), etc. The last step 

corresponds to |so — 3, so — 3) ^ |so — 7/2, so — 7/2). 

Since the individual substeps in Fig. ^(f ) are higher for 
the S^-decreasing transitions (to the right of the plot) 
than for the 5*^ -increasing transitions (to the left), the 
second large-scale step in Fig. ^(e) is higher than the 
first. The reason for this height difference lies in the 
fact that matrix elements for transitions between differ- 
ent spin wave states contain Clebsch-Gordan coefficients: 
Let us consider the two multiplets which give rise to the 
current steps in Fig. ^(f), namely the (sg — 3) and the 
(so — 7/2)-niultiplet. The energy levels of these multi- 
plets are schematically depicted in Fig. ||. The first tran- 
sition between the multiplets which becomes possible as 
the applied voltage is increased is a S*^ -increasing transi- 
tion, namely the |so — 3, so — 4) — > |so — 7/2, so — 7/2)- 
transition at the bottom of the ladder. As the voltage 
is increased further more and more S*^ -increasing tran- 
sitions become possible. (In Fig. ^ they are numbered 
as 1,2,- • ■ ,2so — 6.) The last S'^-increasing transition, 
namely |so - 3, -(sq - 3)) |so - 7/2, -(so - 7/2)), 
lies at the top of the ladder. As we increase the 
applied voltage further the first 5^-decreasing tran- 
sition (|so-3,-(so-4)) ^ |so-7/2,-(so-7/2))) 
becomes possible (at the top of the ladder). As 
we increase the voltage further more and more S^- 
decreasing transitions become possible, the last transi- 
tion (|so — 3, So — 3) —> |so — 7/2, sq — 7/2)) lying at the 
bottom of the ladder. Let us now compare Clebsch- 
Gordan coefficients involved in the different transitions. 
In Fig. H their magnitude is schematically indicated 
by the width of the arrows marking the transitions. 
For example, the Clebsch-Gordan coefficient involved in 
the matrix element for the first S'^-increasing transi- 
tion (|so — 3, So — 4) |so — 7/2, so — 7/2)) as well as 
the first 5^-decreasing transition (|so — 3, — (so — 4)) — > 
|so - 7/2, -(so - 7/2))) in Fig.|(f) (and Fig.|) is of order 
0(l/s), hence these transitions carry very little weight. 
By increasing the applied voltage, so that more tran- 
sitions become accessible, the Clebsch-Gordan coeffi- 
cients for S'^-increasing transitions increase from 0(l/s) 
to 0(1) as we go up the ladder and those for 5*^- 
decreasing transitions likewise increase from 0(l/s) to 
0(1) as we subsequently go down the ladder. Now,, the 
occupation probabilities P, in general, are largeiH for 
lower lying spin wave states than for higher- lying ones. 
Since transitions out of these more strongly-populated 
lower lying states of the spin-(so — 3)-multiplet have a 
Clebsch-Gordan coefficient of order 0(1) if they are S^- 
decreasing, and of order 0(l/s) if they are S^-increasing, 
we conclude that the total weight of all S'^-decreasing 
transitions is larger than that of all iS^-increasing tran- 
sitions. The net result is that the S'^-decreasing step in 
Fig.^(e) is higher than the S'^-increasing step. 

The relative heights of the two large-scale steps in 
Fig. |^(e) are analyzed in Fig. |^, or more precisely, the 



11 



|s„-7/2,-(s„-7/2) > 

IV3,-(s„-3) >^ 2s„'6--*^ 

^> | s„-7/2.-(^-9/2) > 
|s„-3,-(s„-4)>^ ^;>.^ 

> | s„-7/2.-(^-11/2 )> 
|So-3,-(s„-5)>^g'!-8^^ 



|s„-7/2,s„-11/2> 
.n"* |s„-7/2,s„-9/2> 



|s„-7/2,s„-7/2> 



FIG. 8: Transitions from the (so - 3) to the (so - 7/2)- 
multiplet. Arrwos indicate all transitions allowed by the se- 
lection rule that 5^ can change only by ±1/2. The width of an 
arrow schematically indicates the size of the Clebsch-Gordan 
coefficient involved in the matrix element of the correspond- 
ing transition. The energy separations between the various 
levels are not drawn to scale. The order in which transitions 
become possible as the applied voltage is increased is given by 
the numbers in the arrows, which range from 1 to 2so — 6 for 
S^-increasing transitions, and then from 2so — 5 to 4so — 12 
for 5*^ -decreasing transitions. 



height of the current step due to all S'^-decreasing tran- 
sitions, A/decreasing, relative to the total height of the 
current step A/. Fig. |^ confirms that -decreasing tran- 
sitions typically carry significantly more weight than - 
increasing transitions. Moreover, with increasing spin so 
the height of A/docrcasing increases strongly relative to 
A/, so that for large sq no Zeeman splitting is expected. 

Let us now summarize the consequences of the 
nonequilibrium scenario discussed above for the tunnel- 
ing spectra measured by DGR. (i) Firstly, the temper- 
ature in DGR's experiments, namely T « 80 mK, is 
too large for the fine current steps of Fig. ^(c) and (f) 
due to spin wave excitations to have been observable. 
Instead, only the large-scale current steps of Fig. ||(a) 
and (d) would be observable. The observed resonance 
spacing of SE^os ~ 0.2meV [cf. (PI)] indeed does agree 
with that expected for minorit y- ele ctron charging tran- 
sitions and So 1000 [cf. Eq. (^)]. (ii) Secondly, the 
nonequilibrium scenario discussed above can also account 
for the fact (P3) that the vast majority of the observed 
transitions within a given sample shift in energy with a 
similar slope for large magnetic fields, since for sp ^ 1: 
S'^-decreasing transitions carry far more weight than S*^- 
increasing transitions [Fig. m. 



D. Spin and single particle excitations 
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FIG. 9: Illustration of single particle excitations: Exam- 
ples of excited single particle states, which can be reached 
from the ground state |s,s)^ by (a) a majority electron 
or (b) a minority electron entering the grain, or from a 
\s — 1/2, s — 1/2) ^ -state by (c) a majority electron or (d) 
a minority electron leaving the grain. 



states (spin wave states) of the grain lead to resonances 
in tunneling spectra spaced by dmin ~ U /2, which for 
So = 1000 gives a value of « 0.2meV, as observed in ex- 
periments. However, these scenarios are not always able 
to explain the large number of resonances observed, since 
njrj'g given in Eq. (|2^) depends strongly on the threshold 
charging energy and can become as small as two when 
^tiircsh order or smaller than dmin- 

We shall now argue below that the abundance of res- 
onances measured by DGR can be explained by taking 
the analysis one step further, namely by including single- 
particle excitations in addition to spin excitations. For 
the spin excitations we shall henceforth restrict our con- 
siderations to transitions between spin ground states as 
in Sec VB, since we saw in Sec. VC that DGR's ex- 



We saw in Sees. VB (and VC) that a nonequilib- 
rium scenario involving transitions between different spin 



perimental temperature was too high to resolve the fine 
structure due-to spin wave transitions. Furthermore, we 
shall assumeCJ r**' ^ F*"", and hence shall take into ac- 
count only excited single particle states involving a single 
particle- hole pair, i.e. states which can be reached from 
the corresponding ground states by a single tunneling 
transition, namely a majority/minority electron enter- 
ing the iV-electron grain or a majority/minority electron 
leaving the iV -f 1-electron grainHj Fig. ^ illustrates some 
examples. 

We solved Ec|^(|To|) numerically for the parameters of 
Fig. 1 in Ref. and calculated the current; the result 
is shown in Fig.^(b). The current shows 16 steps in the 
plotted region. Note that the calculation of the current 
was done for the same parameters as in Fig. ^(a) where 
only spin excitations were taken into account. The ar- 
rows in (b) mark all additional resonances with respect to 
(a) which arise due to a combination of spin- and single 
particle excitations. Most of the resonances in Fig. ^(b) 
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age via Eq. E3l Fig. 10 shows clearly that the number 
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10: Resonances as function of iJ*J"<=sh and eVB^ for the 



FIG 

same parameters as Fig, 
resonances appear. 



As is increased additional 



are due to tunneling of minority electrons, namely 14 res- 
onances out of 16. This results from the very different 
density of states of minority and majority electrons at 
the Fermi energy [cf. (P4)]. 

Thus, when single-particle excitations are considered 
in addition to nonequilibrium spin accumulation, addi- 
tional resonances appear at higher voltages, so that the 
number of resonances increases significantly. The reso- 
nances are no longer equally spaced, as was the case for 
pure spin excitations, but the average spacing is of the 
same order of magnitude. 

We also investigated the gate-voltage dependence of 
these resonances. Fig. |l0| shows resonances as function 
of bias voltage {cVBl) and threshold charging energy 
^tiucsh jv^Q^g ^Yiat E^^'"^'^ depends linearly on gate voh- 



of resonances increases with increasing E^'^"^^ [cf.(P5)]. 
This is in qualitative agreement with experimental data 
[Fig. 3(a) in Ref. |^, although a detailed understanding 
of the observed gate voltage dependence of resonances 
would require a more systematic examination. 



VI. CONCLUSIONS 

In summary, we have shown that nonequilibrium spin 
and single-particle excitations within our model are able 
to explain most of the experimental data by DGR on tun- 
neling spectroscopy of ultrasmall ferromagnetic grains. 
Nonequilibrium spin excitations lead to a resonance spac- 
ings and an absence of Zeeman splitting as observed in 
experiments. By taking single-particle-excitations into 
account in addition to spin excitations, it is possible to 
satisfactorily explain the large number of resonances ob- 
served in experiment. 
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